Radiative corrections to the Dalitz plot of decays: Contribution of the four-body region 
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We calculate the radiative corrections to the Dalitz plot of Kj^ decays to order (a/ir)(q/Mi), where q is the 
momentum transfer and Mi is the mass of the kaon. We restrict the analysis to the so-called four-body region, 
which arises when no discrimination of real photons is made either kinematically or experimentally. We present 
our results in two ways: the first one with the triple integration over the photon kinematical variables ready to be 
performed numerically and the second one in a fully analytical form. Our results can be useful in experimental 
analyses of the Dalitz plot, by evaluating the model-independent coefficients of the quadratic products of the 
form factors; we provide some numbers as a case example. We find a small, albeit non-negligible, contribution 
from the four-body region to the radiative correction to the total decay rate of decays. 

PACS numbers: 14.40.Df, 13.20.Eb, 13.40.Ks 



There are several inherent difficulties in the analysis of 
radiative corrections (RC) in decays. The absence of first 
principles to evaluate them introduces a model dependence so 
experimental analyses which use them also become model- 
dependent. RC depend on the process characteristics, such 
as the charge assignment of the participating mesons, the 
size of the momentum transfer, and whether real photons can 
be experimentally discriminated or not. RC also depend on 
the observable to be measured. All this requires RC to be 
recalculated everytime the process characteristics and/or the 
observable are changed. Finally, there are difficulties of a 
practical nature: It turns out that the final expressions of RC 
calculations are rather inefficient to use or are long and tedious 
to the point that their use becomes unreliable. 

In a recent paper |Q]] we overcome the above difficulties 
rather satisfactorily. We obtained a model-independent ex- 
pression for the Dalitz plot of Kj^ decays including RC of 
order (a/V)(g/Mi), where q is the momentum transfer and 
Mi is the mass of the kaon. We studied the so-called three- 
body region (TBR) of the Dalitz plot. We thus assumed 
that not only real photons are not detected but also that 
events whose energies do not satisfy the three-body energy- 
momentum conservation restrictions are rejected. 

When real photons cannot be discriminated, however, the 
TBR of the Dalitz plot should be extended to the four-body 
region (FBR). The purpose of the present paper is to extend, 
on the same footing, the calculations of Ref. QJJ] to evaluate the 
four-body contribution of the RC to the Dalitz plot. Although 
this problem keeps a close similarity with the one discussed 
in Ref. [1], it is not possible to find a single rule which 
allows us to adapt one calculation into the other. Thus, we 
need to take a few steps backwards and start the calculation 
at the bremsstrahlung transition amplitude level, namely, by 
squaring such amplitude and then performing the integrals 
over the kinematical variables of the photon restricted to the 



FBR. In order to save a substantial amount of effort we will 
follow the approach implemented in the FBR analysis of RC 
in baryon semileptonic decays ||2|-|4j]. 

For definiteness, let us consider the four-body decay 

K+( Pl ) -> ir°( P2 ) + £+(l) + Vl ( Vv ) + j(k), (1) 

where K + and 7r° denote a positively charged kaon and 
a neutral pion, respectively, £ + and vi denote a positively 
charged lepton (£ + — e + or and its neutrino, respectively, 
and 7 represents a real photon. Here p\ = (E\, pi), p 2 — 
(E 2 ,p 2 ), I = (E,l), p v = (E v ,p v ), and k = (w,k) are the 
four-momenta of K + , 7r°, £ + , v, and 7, respectively. Mi, M 2 , 
and m are the nonzero masses of the first three particles. In 
the center-of-mass frame of K + , Mi = E + E 2 + E v + u) 
and = p 2 + 1 + p, y + k, so E v = E° - uj and p^ = p" - k, 
where E® is the energy and p^ is the three-momentum of the 
neutrino in the nonradiative process. In addition, p 2 , 1, p v , and 
k will also denote the magnitudes of the corresponding three- 
momenta when the expressions involved are not manifestly 
covariant. All other conventions and notation are given in 
Ref. fl]. 

The calculation of bremsstrahlung RC in the FBR simplifies 
because the events in this region have the same amplitude 
M B as in the TBR, given by Eq. (41) of Ref. fT]. The 
bremsstrahlung differential decay rate in the FBR is then 
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In order to perform the integrals over the kinematical variables 
in Eq. (|2), we orient the coordinate axes in such a way that 
£ + is emitted along the +z axis and ir° is emitted in the first 
or fourth quadrant of the (x, z) plane. Therefore, there are 
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five relevant variables of the final state. Two of them are the 
energies E and E 2 , whose allowed values in the FBR are 
given by Eq. (5) of Ref. fl]. The other three variables can 
be grouped into either (k,cos9k,ipk) or (cos 62, cosOk, <fik), 
where k, 9 k , and tp k are the magnitude of the three-momentum 
and the polar and azimuthal angles of the photon, respectively, 
and 62 is the polar angle of ir° . If we define p 2 1 = cos 9 2 = y 
and 1 ■ k = cos 9^ = x, the photon energy is given by [ 1 ] 



where 



F 

2D' 



(3) 



The organization of Eq. (0 allows one to contrast it with 
its counterpart in the TBR, constituted by Eqs. (68)-(73) of 
Ref. where the factors dQ and Ao are defined. The 
counterpart of the first summand on the right-hand side of 
Eq. © is the product AoIq, where the function In, defined 
in Eq. (65) of this reference, contains the infrared divergence. 
Now, Iq is replaced by the infrared-convergent function /of 
defined as 
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F = 2p 2 l{y -y), 
D = E° + Ix + p 2 • k, 
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where 



with 
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By following a standard procedure, the differential decay 
rate in the FBR can be written as 



dr F B BR 



~dn{A I 0F (E,E 2 ) 

7T 



-'of =4 [ — arctanh/3 — 1 



(8) 



P2I 

47T 



and (3 = l/E. All other contributions of |Ms| 2 can be cast 
into a single term |M' B | 2 , whose form is not needed here. Let 
us notice that the variable y falls within the intervals [—1, y ] 
and [—1,1] in the TBR and FBR, respectively. 



J dx J dy J g?<^/c|M' b | 2 >. (6) Performing the integrals in Eq. (O yields 
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where, for simplicity, hereafter we will use the shorthand 
notation f± = f±(q 2 ) unless explicitly noted otherwise. Let 
us stress that there is a one-to-one correspondence between 
the five terms involved in the second expression on the right- 
hand side of Eq. (O and Eqs. (69)-(73) of Ref. yl, so we have 
kept the parallelism between both calculations all along. The 
coefficients C, read 
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The functions A^f (k = 1, ...,11) are the result of 
performing the triple integrals contained in Eq. ©. They have 
the very same structure as their TBR counterparts listed in 
Appendix B of Ref. [1], except for the fact the upper limit of 
integration over the variable y is replaced by 1 . 



Equation (O constitutes our first result: an expression 
for the bremsstrahlung RC to the if,g Dalitz plot to order 
(a/Tr)(q/Mi), restricted to the FBR. At this stage of the cal- 
culation, <ir™ R contains triple integrals over the kinematical 
variables of the photon, which can be performed numerically. 

We can proceed further, however, and provide a fully 
integrated expression that could be more useful in a Monte 
Carlo simulation. For this task, we use previous results 
obtained in the analysis of the FBR of the Dalitz plot of baryon 
semileptonic decays ||2|-|4l]. Accordingly, the analytical form 
of Eq. (|9]l can be expressed as 



dr F B BR 



where 



a 



^drt [^ifI/+I 2 + a 2 fRc(/ + /*) + a 3 f|/-| 2 ], 

(ii) 



An, = A^ Jqf + AlP , i= 1,2,3. (12) 
The functions A^ are defined in Eqs. (17)-(19) of Ref. fl], 
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whereas the functions A\ F ' read 

g 

A 2F ] = -^2(-A 2F -2A3 F -A 5F -2A 6F + A 7F + A 9F 
+ 2A 11F ), (14) 



and 

g 

= ~Ml ^ 3F + AeF " A§F ~ Al0F + AllF ^' (15) 

where the functions Afe F (k — 1, ...,11) correspond to 
the analytical versions of their counterparts A^ of the TBR, 
Eqs. (74)-(84) of Ref. [1]. After some rearrangements, one 
finds in a close analogy with Eqs. (87)-(89) of this reference 
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Equation ( fTTT i constitutes our second result. It is an analyt- 
ical expression for the Dalitz plot of Kj^ decays, restricted 
to the FBR region. It includes RC to order (a/ir)(q/Mi) 
and is model independent. The complete RC to the Dalitz 
plot of K£ decays without the restriction of eliminating real 
photons are given by adding <fT™ R to <fT™ R given by Eq. (93) 
of Ref. 01. 



An application of our formulas in a Monte Carlo analysis 
requires the numerical evaluation of the functions Ai and 
Aip all over the Dalitz plot. This can be done by tracing 
a lattice in the allowed kinematical region and constructing 
arrays with these functions evaluated at given points (E, E 2 ). 
The arrays should feed the Monte Carlo simulation as a matrix 
multiplication, with the form factors and the energies (E, E 2 ) 
varied in each step of the simulation. In Tables U and [TT] 
we display the numerical values of (a/n)Ai and (a/Tv)A±F, 
the latter in boldface characters, for K^ 3 and KT 3 decays, 
respectively, as a case example. 



As a further application, we can evaluate the total decay 
rate of Kf 3 including RC from both the TBR and FBR. It can 



be written as 
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where I(X+) and I F (A 

the energies E and E 2 for the TBR and FBR 
Their expressions read 
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) involve the double integration over 
respectively. 
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TABLE I: Radiative corrections {a/it)A\ x 10 and (a/ir)AiF x 10 for K + — s> it + e + + v & decay. The entries corresponding to the latter 
are marked in boldface characters. E and Ei are given in GeV. 



E 2 \E 


0.0123 


0.0370 


0.0617 


0.0864 


0.1111 


0.1358 


0.1604 


0.1851 


0.2098 


0.2592 


0.1533 


0.1880 


0.1462 


0.0668 


-0.0286 


-0.1220 


-0.1949 


-0.2246 


-0.1726 


0.2468 


0.1028 


0.1810 


0.1580 


0.0902 


0.0011 


-0.0905 


-0.1654 


-0.1998 


-0.1537 


0.2345 


0.0701 


0.1578 


0.1522 


0.0962 


0.0150 


-0.0718 


-0.1445 


-0.1792 


-0.1354 


0.2222 


0.0490 


0.0948 


0.1429 


0.0989 


0.0261 


-0.0551 


-0.1249 


-0.1590 


-0.1168 


0.2098 


0.0339 


0.0599 


0.1321 


0.1004 


0.0363 


-0.0391 


-0.1056 


-0.1389 


-0.0977 


0.1975 


0.0228 


0.0383 


0.0662 


0.1014 


0.0461 


-0.0233 


-0.0863 


-0.1186 


-0.0779 


0.1851 


0.0145 


0.0236 


0.0375 


0.0743 


0.0558 


-0.0075 


-0.0670 


-0.0979 


-0.0568 


0.1728 


0.0084 


0.0134 


0.0203 


0.0343 


0.0654 


0.0083 


-0.0474 


-0.0769 


-0.0333 


0.1604 


0.0042 


0.0065 


0.0095 


0.0149 


0.0281 


0.0243 


-0.0275 


-0.0550 


-0.0019 


0.1481 


0.0014 


0.0021 


0.0030 


0.0045 


0.0076 


0.0170 


-0.0070 


-0.0316 





TABLE II: Radiative corrections (a/-iv)Ai x 100 and (a/ir)Aip x 100 for K + — s> 7r° + fi + + decay. The entries corresponding to the 
latter are marked in boldface characters. E and E2 are given in GeV. 



E 2 \E 


0.1131 


0.1280 


0.1429 


0.1578 


0.1727 


0.1876 


0.2025 


0.2174 


0.2322 


0.2480 










-0.3664 


-0.2633 


-0.2285 


-0.1879 


-0.1003 


0.2361 




-0.0168 


-0.0295 


-0.0624 


-0.1026 


-0.1406 


-0.1649 


-0.1584 


-0.0861 


0.2242 


0.0397 


0.0356 


0.0096 


-0.0301 


-0.0760 


-0.1188 


-0.1472 


-0.1437 


-0.0725 


0.2123 


0.0345 


0.0355 


0.0139 


-0.0224 


-0.0657 


-0.1068 


-0.1341 


-0.1303 


-0.0582 


0.2004 


0.0148 


0.0237 


0.0091 


-0.0214 


-0.0600 


-0.0975 


-0.1222 


-0.1170 


-0.0438 


0.1885 


0.0238 


0.0072 


0.0006 


-0.0231 


-0.0560 


-0.0890 


-0.1106 


-0.1036 




0.1766 


0.0127 


0.0651 


-0.0096 


-0.0260 


-0.0527 


-0.0809 


-0.0991 


-0.0899 




0.1647 


0.0067 


0.0213 


-0.0207 


-0.0294 


-0.0498 


-0.0729 


-0.0874 


-0.0761 




0.1528 


0.0030 


0.0087 


0.0213 


-0.0332 


-0.0470 


-0.0649 


-0.0758 


-0.0652 




0.1409 


0.0008 


0.0024 


0.0052 


0.0117 


0.0498 


-0.0574 


-0.0663 
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where f+ + ™ (q 2 ) has been expanded linearly in q 2 , with the 

slope parameter A+ = 0.0328 ± 0.0033 10. 

By using more refined numerical integration routines than 



in Ref. [ 1], we find that the uncorrected coefficients are h y 



(0) 



0.0965, h { ° ] = 0.3568, and h 2 u > = 0.5279. The inclusion of 
RC yields h = 0.0958, hi = 0.3532, and h 2 = 0.5209 
from the TBR and h 0F = 0.0005, h lF = 0.0031, and h 2F = 
0.0062 from the FBR. With the value of A + quoted above, we 



,(0) 



find that RC from the TBR induce a decrease of 0.8% in the 
decay rate and the combined TBR and FBR effect produces 
an overall decrease of 0.26%, which is an important effect. 

In summary, we have performed a detailed analysis of the 
computation of RC from the FBR of the Dalitz plot to order 
(a/ir)(q/Mi), where q is the momentum transfer and Mi 
is the mass of the kaon. At this order of approximation, 
the expression obtained is model independent; it has been 
organized in such a way that its use in a Monte Carlo 
simulation is straightforward. The effects of the FBR are 
quite important when no discrimination of real photons can be 
made either kinematically or by direct detection; they produce 
a perceptible decrease in the total decay rate. Our results may 
constitute an important theoretical input in the determination 
of V us from K13 decays. 
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